Abstract
The network topology in disordered materials is an important structural descriptor for understanding the nature of disorder that is usually hidden in pairwise correlations. Here, we compare the covalent network topology of liquid and solidified silicon (Si) with that of silica (SiO2) on the basis of the analyses of the ring size and cavity distributions and tetrahedral order. We discover that the ring size distributions in amorphous (a)-Si are narrower and the cavity volume ratio is smaller than those in a-SiO2, which is a signature of poor amorphous-forming ability in a-Si. Moreover, a significant difference is found between the liquid topology of Si and that of SiO2. These topological features, which are reflected in diffraction patterns, explain why silica is an amorphous former, whereas it is impossible to prepare bulk a-Si. We conclude that the tetrahedral corner-sharing network of AX2, in which A is a fourfold cation and X is a twofold anion, as indicated by the first sharp diffraction peak, is an important motif for the amorphous-forming ability that can rule out a-Si as an amorphous former. This concept is consistent with the fact that an elemental material cannot form a bulk amorphous phase using melt quenching technique.
Similar content being viewed by others
Introduction
The absence of translational periodicity and symmetry, and the rich structural complexity make it difficult to understand the order within disorder1,2 in disordered materials. The advent of advanced instrumentation and measurement protocols makes it feasible to use quantum beam diffraction (X-ray diffraction (XRD) and neutron diffraction (ND)) techniques to reveal the structure of disordered materials at synchrotron and/or neutron sources3,4,5,6.
Amorphous (a)-silicon (Si) and silica (SiO2) are the most typical and important disordered materials in both fundamental and technological research studies. In particular, SiO2 is a canonical amorphous-former, whereas it is possible to synthesize a-Si only in a thin film owing to its poor amorphous-forming ability. The short-range structural unit of these amorphous materials is a tetrahedron, SiSi4 in a-Si and SiO4 in a-SiO2, and the formation of their networks is governed by the corner-sharing of tetrahedra. This corner-sharing motif is within Zachariasen’s classification7 for the glass formation of oxide materials, but elemental amorphous materials do not follow the rule, because it is impossible to obtain bulk amorphous silicon as mentioned above.
The structures of these materials have been studied by diffraction techniques. Laaziri et al. reported high-quality X-ray diffraction data for a-Si with a high-real space resolution for precisely determining the coordination number8. The structure of a-SiO2 has been widely studied by both X-ray9,10 and neutron diffraction4,11,12. The most important feature in the diffraction data of a-SiO2 is that the first sharp diffraction peak (FSDP)11,12,13,14 is observed in both X-ray and neutron diffraction data, whereas the second diffraction peak, the so-called principal peak (PP)12,14, can be observed in only the neutron diffraction data because this peak reflects the packing of oxygen atoms12,15, which is sensitive to neutrons. The origin of the FSDP has been discussed for a long time. The FSDP was first discussed in 197616, although it seems that the name “FSDP” was first used by Phillips in 198117. The interpretation of diffraction peaks including the FSDP was attempted in the 1980s13,17, as discussed in details in several papers18,19,20,21,22. It is known that the FSDP of a-SiO2 is related to the formation of the random network model proposed by Zachariasen7, which was extended to silicate glasses by Greaves23 and recently revised by Mei et al., as illustrated in Ref. 21. It was demonstrated that intermediate-range ordering arises from the periodicity of boundaries between successive cages in the network formed by the connection of regular SiO4 tetrahedra with shared oxygen atoms at the corners associated with the formation of a ring structure and a large cavity21,22. The second maximum, PP, reflects the size of the local-network-forming motif, whereas the FSDP indicates the arrangement of these motifs in an intermediate range according to Zeidler and Salmon24. Another interpretation of the FSDP has recently been proposed by Shi and Tanaka, who discussed local tetrahedral ordering in covalent liquids and glasses25 and they concluded that a-Si has an FSDP. Moreover, an FSDP was found in metallic glasses26, although Price et al. implied that some diffraction peaks from amorphous alloys are not FSDPs13.
The investigation of the behavior of the FSDP and PP under high pressures is important to understand the nature of intermediate-range ordering in disordered materials. Figure 1a shows in situ neutron structure factors, S(Q), for a-SiO2 under high pressures reported by Zeidler et al.27. It was found that the FSDP diminishes with increasing pressure, simultaneously with a peak shift to a higher Q. On the other hand, the PP becomes sharp with increasing pressure, suggesting that the oxygen packing fraction increases12,15 with the decrease of cavity volume under high pressures. This is an important benchmark for the modification of the intermediate-range ordering of a-SiO2 under cold compression. Onodera et al. recently reported the unusual behaviour of the FSDP after hot compression2; in the X-ray S(Q) shown in Fig. 1b, they have observed the evolution of FSDP at 7.7 GPa at a temperature higher than 400 °C. These diffraction data suggest that FSDP is very sensitive to pressure and temperature, while the PP position is insensitive to the density change.
In this article, we apply several topological techniques to analyze the ring size and cavity distributions, and tetrahedral order of crystalline, amorphous, and liquid Si and SiO2 to reveal the network topology2,28,29,30 for understanding the diffraction peaks in disordered materials with a special focus on the FSDP and PP.
Methods
Structure modeling
Atomistic models of liquid (l)-SiO2 and a-SiO2 were obtained by classical molecular dynamics (MD) simulation and MD–reverse Monte Carlo (RMC) modeling29, respectively. The atomistic models for l- and a-SiO2 were generated by combined classical MD simulation–RMC modeling. The MD model for l-SiO2 was adopted from the literature31. The MD simulation for a-SiO2 was performed using the large-scale atomic/molecular massively parallel simulator (LAMMPS) package32 within the NVT ensemble. The interactions were described by pair potentials with short-range Born–Mayer repulsive and long-range Coulomb terms, i.e.,
where r is the distance between atoms i and j, Bij and ρij define the magnitude and softness of the Born–Mayer terms, respectively, Zi is the effective charge on atom i (ZSi = 2.4, ZO = -1.2), e is the elementary charge, and ε0 is the permittivity of vacuum. The Bij values were 21.39 × 10−16 J (Si–O), 0.6246 × 10−16 J (O–O) or zero (Si–Si); the ρij values were 0.174 Å (Si–O), 0.362 Å (O–O) or zero (Si–Si). As initial configuration, 3000 (Si, 1000; O, 2000) atoms were randomly distributed in a cubic cell with a side length of 35.66 Å. The cell had a number density of 0.06615 Å-3. The simulation used periodic boundary conditions, and the long-range Coulomb interactions were treated by using the Ewald summation. A time step of 1 fs was used in the Verlet algorithm. The simulation temperature was maintained at 4000 K for 20,000 time steps, then the temperature was reduced to 300 K over 200,000 time steps. Finally, the system was equilibrated at 300 K for 50,000 time steps. A Nosé-Hoover thermostat was employed to control the temperature. After MD simulations, the configurations obtained for l- and a-SiO2 were refined by RMC modeling with constraints on the coordination numbers and on the O–Si–O bond angle distribution, in order to prevent the formation of an unfavorable disordered structure.
Atomistic models of l-Si and a-Si were obtained by RMC modeling29 and combined classical MD simulation–RMC modeling, respectively.
The RMC model for l-Si (1770 K, 5000 particles) was obtained by the RMC + + code33 and based on the X-ray S(Q) for l-Si29. The number density was 0.055 Å−3, that is consistent with a bulk mass density of 2.57 g cm−334.
The model for a-Si was obtained by MD simulation, followed by RMC refinement. The MD simulation was performed using the LAMMPS code. The modified Tersoff potential 35 based on the three-body Tersoff potential 36,37 was used for describing interatomic interactions. In the modified Tersoff potential, the total energy E is written as follows:
where E is decomposed into bond energies Vij between atoms i and j. rij is the distance between atoms i and j, θijk is the angle confined by the bonds between ij and ik. The function Aexp(-λ1rij) and bijBexp(− λ2rij) represent a repulsive and an attractive term, respectively. The extra term fC is merely a smooth cutoff function, to limit the range of the potential. bij represents a measure of the bond order. The term 1 + ξij corresponds to the coordination number of atom i, where ij bond and other bonds are counted as 1 and ξij respectively. The term gmod(θ) is the modified angular-dependent term. The potential parameters, A, B, λ1, λ2, n, α, m, c1, c2, c3, c4, c5, h, R1, and R2, are given in Table 1. The simulation box was a cubic cell with a side length of 50.00 Å. In the MD simulation, 6256 Si atoms were placed in the box as an initial configuration with in the NVT ensemble. A time step of 1 fs was used in the Verlet algorithm. The atomic configuration was initialized at random and the system was equilibrated at 3000 K for 500,000 steps. Then it was cooled to 300 K during 5,000,000 steps and annealed at 300 K for 500,000 steps. After the MD simulation, the configuration obtained was refined by RMC simulation with constraints on the coordination number and the Si–Si-Si bond angle distribution, in order to avoid formation of unfavorable disordered structures. The model reproduces the experimental X-ray S(Q).
For l-Si, additional principles (FP)MD simulations within the framework of density functional theory were performed for a 64-atom Si supercell with periodic boundary conditions. The calculations were performed using the projector-augmented wave method38 and the generalized-gradient approximation with the exchange correlation functional of Perdew, Burke, and Ernzerhof39. The electronic wave functions were expanded in a plane-wave basis with an energy cutoff of 400 eV at the Γ point in the Brillouin zone. Atomic configurations selected from the FPMD trajectory generated using the Vienna Ab initio Simulation Package40 were used to validate the atomistic model of l-Si.
Ring size distribution analysis
The ring size distribution was calculated using the R.I.N.G.S. code41,42 on the basis of the King43 and primitive44,45 criteria. The first coordination distances were set to 2.5 Å and 3.0 Å for a-Si and l-Si, and to 1.9 Å and 2.4 Å for a-SiO2 and l-SiO2, respectively.
Cavity volume calculation
Cavity volume analyses were performed using the pyMolDyn code46. The code can calculate three different types of cavity, domain, center-based (Voronoi), and surface-based cavities. We calculated surface cavity volumes with a cut off distance rc of 2.5 Å.
Tetrahedral order parameter q
The tetrahedral order parameter for the Si-centered hyper-tetrahedra and SiO4 tetrahedra are defined as47
where θijk is the angle formed between the central Si atom j and its neighboring Si or O atoms i and k. This parameter was designed to give a value of unity for a regular tetrahedron and a mean value of zero for a perfect gas.
Results and discussion
The X-ray structure factor, S(Q), for a-Si8 is shown in Fig. 2a together with that of l-Si (1770 K)29. We observe Q2 (PP) and Q3 at QrA-X ~ 5 and 8.5, respectively, and no Q1(FSDP) is observed. Note that scattering vector Q is scaled by multiplying by rA-X (first coordination distance)6,11,12,13,14,15,18,22,24,29,48 obtained by a Fourier transform of S(Q). S(Q) for l-Si is different from that for a-Si owing to the increased coordination number, NSi-Si, from approximately 3.9 (amorphous)8 to 5.7 (liquid)29 associated with the significant density increase from 2.30 g cm−3 (amorphous) to 2.57 g cm−3 (liquid). This behavior was also confirmed in previous FPMD studies49,50 and is consistent with the fact that a-Si is a semiconductor and l-Si is a metallic liquid. X-ray29,30 and neutron30 structure factors, S(Q), for a-SiO2 are shown in Fig. 2b. As mentioned in the previous section, Q1 (FSDP) is observed for a-SiO2 in both the X-ray and neutron S(Q), but Q2 (PP) is visible only in the neutron S(Q), because it reflects the packing fraction of oxygen atoms12,15 since neutrons are sensitive to O–O correlations but X-rays are more sensitive to Si–Si correlations. It is worth mentioning that the X-ray S(Q) of l-SiO2 (2019 K) is identical to that of g-SiO2. Q1 (FSDP) in the X-ray S(Q) shown in Fig. 2b is prominent in l-SiO251, suggesting that Si–O covalent bonds are strong even in the liquid (2019 K). This behavior is consistent with the small differences in density (amorphous: 2.21 g cm-3, liquid: 2.10 g cm-3) and Si–O coordination number, NSi-O (amorphous: 4.0, liquid: 3.9) between them. Considering the melting point of Si (1687 K)52 and SiO2 (1996 K)53, we have noted that the S(Q) of l-Si markedly changes at only 100 K above the melting point, although the temperature-dependent change in the diffraction pattern is small in l-Si32 (see Fig. S1). The marked change in structure above the melting point suggests a large difference in configuration entropy, which is evidence of the poor amorphous-forming ability. On the contrary, the S(Q) of l-SiO2 is similar to that of g-SiO251. We assume that the change in the S(Q) is also an indicator of the amorphous-forming ability.
Recently, Shi and Tanaka have reported a comparison of S(Q) of several disordered materials with a tetrahedral motif25, by scaling rA-X and rA-A (cation–cation distance). We also display the S(Q) of Si and SiO2 in Fig. 3. As can be seen in the figures, a-Si exhibits an FSDP when the S(Q) is scaled by rA-A, but not when scaled by rA-X. The scaling of Q is still an open question, but, recently Salmon and Zeidler have reported that a-Si does not have an FSDP12.
To understand the relationship between the diffraction peaks and the topology, we calculated ring size distributions on the basis of King and primitive criteria. The primitive and King ring size distributions for crystalline (c)-Si and SiO2 (cristobalite) calculated from crystal structures, and for structure models for liquid and amorphous Si and SiO2 obtained by computer simulations, are shown in Fig. 4 and Fig. S2, respectively. Note that the ring criterion does not affect the overall shape of the distributions. The crystalline phases of both Si and SiO2 exhibit only sixfold rings, in which the ring comprises six silicon atoms in c-Si, whereas it comprises six silicon atoms and six oxygen atoms (AXAX rings) in c-SiO2. The ring size distributions of a-Si and a-SiO2 are broad, and the sixfold rings are dominant for both a-Si and a-SiO2. The ring size distributions of a-Si are consistent with the result of a previous study reported by Opletal et al.54. It is stressed that the ring size distributions of l-Si are significantly different from those of a-Si. The large fraction of threefold rings in the liquid is due to the increased coordination number associated with a significant increase in the density of the liquid. This behavior of the ring size distribution is in line with the difference in diffraction data shown in Fig. 2. On the other hand, the ring size distributions of l-SiO2 (2373 K) are identical to those of a-SiO2, which is consistent with the similarity of the density and coordination number between them. The broad ring size distributions of a-SiO2 and l-SiO2 are a signature of topological disorder according to Gupta and Cooper55, and we previously hypothesized that the width of the ring size distribution is an indicator of the amorphous-forming ability28,56 when short-range structures are the same. Indeed, it appears that a-Si exhibits a narrower ring size distribution than a-SiO2, demonstrating that topological order/disorder is a suitable descriptor for understanding the amorphous-forming ability.
The formation of rings is an important structural feature in covalent amorphous materials; hence, it results in the generation of cavities (empty spaces). We visualize a cavity (highlighted in green) of a-SiO2 together with the cavity size distributions in a-SiO2 and l-SiO2 (2373 K) in Fig. 5. The cavity volume ratio of a-SiO2 is approximately 33%29. The large cavity volume ratio is presumed to be a signature of covalent amorphous materials. Note that the cavity volume ratio of l-SiO2 is comparable to that of a-SiO2 owing to the small difference in density between them. The cavity size distributions in SiO2 indicate that both a- and l-SiO2 have a large cavity in disordered structure. Figure 6 visualizes cavities in a-Si and l-Si together with the cavity size distributions. The cavity volume ratio of 22% for a-Si is much smaller than that of 33% for a-SiO2 shown in Fig. 5, but we observe a significant reduction in cavity volume ratio in l-Si (from 22% in a-Si to 2% in l-Si), which is associated with the increase in density (2.30 to 2.57 g cm−3) and coordination number NSi-Si (3.9 to 5.7). It is worth mentioning that the cavity size distributions indicated that a large cavity in a-Si is squeezed and only small cavities are observed in the liquid. This behavior suggests a poor amorphous-forming ability of Si, which is in line with the fact that a-Si can only be formed in a thin film.
Shi and Tanaka have recently focused on the symmetry of SiSi4 tetrahedra in a-Si and SiSi4 hyper-tetr ahedra in a-SiO225. We apply this analysis to crystalline and amorphous Si and SiO2. Figure 7 shows q values of a series of Si and SiO2. Both c-Si and c-SiO2 have perfect tetrahedra as indicated by the average q value of 1. A comparison between a-Si and a-SiO2 demonstrates that SiSi4 hyper-tetrahedra in a-SiO2 are highly distorted in comparison with regular SiSi4 tetrahedra in a-Si, as revealed by the average q values of a-SiO2 being much smaller than that for a-Si of 0.90. In addition, the distribution of q values is broader in a-SiO2. It is also suggested that hyper-tetrahedra in l-SiO2 are slightly distorted in comparison with those in a-SiO2, whereas the average q value of l-Si is the smallest owing to the formation of fivefold and sixfold Si in the liquid. The profile of the q-parameter for l-Si is consistent with the profile reported in a previous FPMD study49,50. Therefore, the structural motif in the liquid phase is completely different between Si and SiO2 as indicated by the behavior of the density, since the density difference between a-SiO2 and l-SiO2 is small, but large between a-Si and l-Si.
Thus, we have revealed the differences in terms of diffraction peaks and the topology by a combination of diffraction measurements and topological analyses. We show that the diffraction data and the topology for a-Si and l-Si are very different. This result is in sharp contrast to SiO2, in which both the diffraction pattern and topology are identical for a-SiO2 and l-SiO2. We are confident that this behavior is the reason why SiO2 is a glass former and Si is not. In other words, the tetrahedral corner-sharing network of AX2, in which A is a fourfold cation and X is a twofold anion, as indicated by the FSDP, is an important motif for the amorphous-forming ability that can rule out a-Si as a good amorphous former. This concept is consistent with the fact that it is impossible to form an elemental bulk amorphous material using melt quenching technique57.
As mentioned in the introduction, the FSDP and PP are suitable descriptors of an amorphous network. The typical behavior of the FSDP and PP under high pressure shown in Fig. 1, in which the FSDP diminishes with a peak shift to a higher Q and the PP becomes sharp with increasing pressure at room temperature, without any peak shifts (Fig. 1a). Onodera et al. have recently reported the evolution of the FSDP with a peak shift in a-SiO2 owing to an LDA-HDA transition induced by heating under high pressure (see Fig. 1b)2. No such significant peak shift can be observed in a-Si under pressures owing to the very strong fully tetrahedral covalent network associated with the absence of an FSDP (note, however, that a first-order-like transformation to high-density a-Si occurs above 10 GPa58,59,60,61). Therefore, it is suggested that the first peak in the S(Q) of a-Si can be assigne d to a PP, because the peak does not show any shifts under high pressures. This interpretation supports that FSDP is a signature of amorphous formin g ability in amorphous network system, although FSDP is not a signature of network29. On the other hand, the FSDP of liquid phosphorus diminishes under high pressures and temperatures, which is associated with the network formation of P4 molecules by transition62,63.
In this article, we demonstrate that topological analyses are powerful tools for understanding the diffraction peaks of disordered materials. Our conclusion fully supports the discussion of Price et al. in 198813. Understanding the diffraction peaks of disordered materials via a series of topological analyses will give rise to the capability to forge a new path for designing novel functional disordered materials.
References
Salmon, P. S. Order within disorder. Nat. Mater. 1, 87–88 (2002).
Onodera, Y. et al. Structure and properties of densified silica glass: Characterizing the order within disorder. NPG Asia Mater. 12, 85 (2020).
Fischer, H. E., Barnes, H. E. & Salmon, P. S. Neutron and X-ray diffraction studies of liquids and glasses. Rep. Prog. Phys. 69, 233–299 (2006).
Greaves, G. N. & Sen, S. Inorganic glasses, glass-forming liquids and amorphizing solids. Adv. Phys. 56, 1–166 (2007).
Kohara, S. & Salmon, P. S. Recent advances in identifying the structure of liquid and glassy oxide and chalcogenide materials under extreme conditions: A joint approach using diffraction and atomistic simulation. Adv. Phys. X 1, 640–660 (2016).
Ohara, K., Onodera, Y., Murakami, M. & Kohara, S. Structure of disordered materials under ambient to extreme conditions revealed by synchrotron x-ray diffraction techniques at SPring-8–Recent instrumentation and synergic collaboration with modelling and topological analyses. J. Phys. Cond. Matter 33, 383001 (2021).
Zachariasen, W. H. The atomic arrangement in glass. J. Am. Chem. Soc. 54, 3841–3851 (1932).
Laaziri, K. et al. High resolution radial distribution function of pure amorphous silicon. Phys. Rev. Lett. 82, 3460–3463 (1999).
Poulsen, H. F. et al. Amorphous silica studied by high energy X-ray diffraction. J. Non-Cryst. Solids 188, 63–74 (1995).
Kohara, S. & Suzuya, K. Intermediate-range order in vitreous SiO2 and GeO2. J. Phys. Condens. Matter 17, 77–86 (2005).
Wright, A. C. Neutron scattering from vitreous silica. V. The structure of vitreous silica: What have we learned from 60 years of diffraction studies?. J. Non-Cryst. Solids 179, 84–115 (1994).
Salmon, P. S. & Zeidler, A. Ordering on different length scales in liquid and amorphous materials. J. Stat. Mech. Theory E 2019, 114006 (2019).
Price, D. L. et al. Intermediate-range order in glasses and liquids. J. Phys. C Solid State 21, L1069–L1072 (1988).
Salmon, P. S., Martin, R. A., Mason, P. E. & Cuello, G. J. Topological versus chemical ordering in network glasses at intermediate and extended length scales. Nature 435, 75–78 (2005).
Salmon, P. S. In Magma under Pressure: Advances in High-Pressure Experiments on Structure and Properties of Melts (eds Kono, Y. & Sanloup, C.) 347 (Elsevier, 2018).
Wright, A. C. & Leadbetter, A. J. Diffraction studies of glass structure. Phys. Chem. Glasses 17, 122–145 (1976).
Phillips, J. C. Topology of covalent non-crystalline solids II: Medium-range order in chalcogenide alloys and a-Si(Ge). J. Non-Cryst. Solids 43, 37–77 (1981).
Elliott, S. R. Medium-range structural order in covalent amorphous solids. Nature 354, 445–452 (1991).
Salmon, P. S. Real space manifestation of the first sharp diffraction peak in the structure factor of liquid and glassy materials. Proc. R. Soc. Lond. Ser. A 445, 351–365 (1994).
Gaskell, P. H. & Wallis, D. J. Medium-range order in silica, the canonical network glass. Phys. Rev. Lett. 76, 66 (1996).
Mei, Q. et al. Intermediate range order in vitreous silica from a partial structure factor analysis. Phys. Rev. B 78, 144204 (2008).
Kohara, S. et al. Atomic and electronic structures of an extremely fragile liquid. Nat. Commun. 5, 5892–5898 (2014).
Greaves, G. N. EXAFS and the structure of glass. J. Non-Cryst. Solids 71, 203–217 (1985).
Zeidler, A. & Salmon, P. S. Pressure-driven transformation of the ordering in amorphous network forming materials. Phys. Rev. B 93, 214204 (2016).
Shi, R. & Tanaka, H. Distinct signature of local tetrahedral ordering in the scattering function of covalent liquids and glasses. Sci. Adv. 5, 3194 (2019).
Ma, D., Stoica, A. D. & Wang, X.-L. Power-law scaling and fractal nature of medium-range order in metallic glasses. Nat. Mater. 8, 30–34 (2009).
Zeidler, A. et al. High-pressure transformation of SiO2 glass from a tetrahedral to an octahedral network: a joint approach using neutron diffraction and molecular dynamics. Phys. Rev. Lett. 113, 135501 (2014).
Akola, J. et al. Network topology for the formation of solvated electrons in binary CaO–Al2O3 composition glasses. Proc. Natl. Acad. Sci. USA 110, 10129–10134 (2013).
Onodera, Y. et al. Understanding diffraction patterns of glassy, liquid and amorphous materials via persistent homology analyses. J. Ceram. Soc. Jpn. 127, 853–863 (2019).
Onodera, Y. et al. Origin of the mixed alkali effect in silicate glass. NPG Asia Mater. 11, 75 (2019).
Takada, A., Richet, P., Catlow, C. R. A. & Price, G. D. Molecular dynamics simulations of vitreous silica structures. J. Non-Cryst. Solids 345&346, 224–229 (2004).
Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).
McGreevy, R. L. & Pusztai, L. Reverse Monte Carlo simulation: A new technique for the determination of disordered structures. Mol. Simul. 1, 359–367 (1988).
Higuchi, K. et al. Precise measurement of density and structure of undercooled molten silicon by using synchrotron radiation combined with electromagnetic levitation technique. Meas. Sci. Technol. 16, 381–385 (2005).
Kumagai, T., Izumi, S., Hara, S. & Sakai, S. Development of bond-order potentials that can reproduce the elastic constants and melting point of silicon for classical molecular dynamics simulation. Comput. Mater. Sci. 39, 457–464 (2007).
Tersoff, J. New empirical approach for the structure and energy of covalent systems. Phys. Rev. B 37, 6991–7000 (1988).
Tersoff, J. Empirical interatomic potential for silicon with improved elastic properties. Phys. Rev. B 38, 9902–9905 (1988).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).
Perdew, J. P., Burke, K. & Ernzerh, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Kresse, G. & Furthmuller, J. Efficiency of ab Initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Roux, S. L. & Jund, P. Ring statistics analysis of topological networks: New approach and application to amorphous GeS2 and SiO2 systems. Comput. Mater. Sci. 49, 70–83 (2010).
Roux, S. L. & Jund, P. Erratum: Ring statistics analysis of topological networks: New approach and application to amorphous GeS2 and SiO2 systems [Comput. Mater. Sci. 49 (2010) 70–83]. Comput. Mater. Sci. 50, 1217 (2011).
King, S. V. Ring configurations in a random network model of vitreous silica. Nature 213, 1112–1113 (1967).
Goetzke, K. & Klein, H. J. Properties and efficient algorithmic determination of different classes of rings in finite and infinite polyhedral networks. J. Non-Cryst. Solids 127, 215–220 (1991).
Yuan, X. & Cormack, A. N. A computational multiscale strategy to the study of amorphous materials. Comput. Mater. 24, 343–360 (2002).
Heimbach, I. et al. pyMolDyn: Identification, structure, and properties of cavities/vacancies in condensed matter and molecules. J. Comput. Chem. 38, 389–394 (2017).
Errington, J. R. & Debenedetti, P. G. Relationship between structural order and the anomalies of liquid water. Nature 409, 318–321 (2001).
Benmore, C. J. et al. Intermediate range chemical ordering in amorphous and liquid water, Si, and Ge. Phys. Rev. B 72, 132201 (2005).
Štich, I., Car, R. & Parrinello, M. Structural, bonding, dynamical, and electronic properties of liquid silicon: An ab initio molecular-dynamics study. Phys. Rev. B 44, 4262–4274 (1991).
Morishita, T. How does tetrahedral structure grow in liquid silicon upon supercooling?. Phys. Rev. Lett. 97, 165502 (2006).
Skinner, L. B. et al. A time resolved high energy X-ray diffraction study of cooling liquid SiO2. Phys. Chem. Chem. Phys. 15, 8566–8572 (2013).
Lide, D. R. (ed.) Handbook of Chemistry and Physics 83rd edn. (CRC Press, 2002).
Wright, A. C. Defect-free vitreous networks: The idealized structure of SiO2 and related glasses. In Defects in SiO2 and Related Dielectrics: Science and Technology 1–35 (Springer, 2000).
Opletal, G., Petersen, T. C., Snook, I. K. & McCulloch, D. G. Modeling of structure and porosity in amorphous silicon systems using Monte Carlo methods. J. Chem. Phys. 126, 214705 (2007).
Gupta, P. K. & Cooper, A. R. Topologically disordered networks of rigid polytopes. J. Non-Cryst. Solids 123, 14–21 (1990).
Kohara, S. et al. Relationship between topological order and glass forming ability in densely packed enstatite and forsterite composition glasses. Proc. Natl. Acad. Sci. USA 108, 14780–14785 (2011).
Hattori, T. et al. Does bulk metallic glass of elemental Zr and Ti exist?. Phys. Rev. Lett. 96, 255504 (2006).
Morishita, T. Structural, electronic, and vibrational properties of high-density amorphous silicon: A first-principles molecular-dynamics study. J. Chem. Phys. 130, 194709 (2009).
Morishita, T. High density amorphous form and polyamorphic transformations of silicon. Phys. Rev. Lett. 93, 055503 (2004).
Daisenberger, D. et al. High-pressure x-ray scattering and computer simulation studies of density-induced polyamorphism in silicon. Phys. Rev. B 75, 224118 (2007).
McMillan, P. F., Wilson, M., Daisenberger, D. & Machon, D. A density-driven phase transition between semiconducting and metallic polyamorphs of silicon. Nat. Mater. 4, 680–684 (2005).
Katayama, Y. et al. A first-order liquid–liquid phase transition in phosphorus. Nature 403, 170–173 (2000).
Morishita, T. Liquid-liquid phase transitions of phosphorus via constant-pressure first-principles molecular dynamics simulations. Phys. Rev. Lett. 87, 105701 (2001).
Acknowledgements
This research was supported by the TIA collabo rative research program "Kakehashi", Grant Number TK19-004 (to S.K. and H.M.), and JSPS KAKENHI Grant Numbers 20H05878 (to M.S., S.K., and K.H.), 20H05881 (to S.K., Y.O., A.H., and K.K.), 20H05884 (to M.S.), 20H04241 (to S.K., M.S., Y.O., and A.H.), 19K05648 (to Y.O.), and 19K22072 (to S.K., Y.O, and H.M.). Discussion with László Pusztai is gratefully appreciated.
Author information
Authors and Affiliations
Contributions
S.K. conceived and designed the research. S.K., M.S., T.M., and Y.O. performed the simulations and modeling. M.S., S.K., H.M., A.H., T.M., K.K., and K.H. analyzed the results. The manuscript was written by S.K. with the input of all the authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kohara, S., Shiga, M., Onodera, Y. et al. Relationship between diffraction peak, network topology, and amorphous-forming ability in silicon and silica. Sci Rep 11, 22180 (2021). https://doi.org/10.1038/s41598-021-00965-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-00965-5
This article is cited by
-
Machine learning molecular dynamics reveals the structural origin of the first sharp diffraction peak in high-density silica glasses
Scientific Reports (2023)
-
Siliceous zeolite-derived topology of amorphous silica
Communications Chemistry (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.